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Abstract 

In this paper we propose a mechanism of natural preheating of our universe by introducing 
an inflaton field dependent mass term for the gravitational wave for a specific class of massive 
gravity theory. For any single field inflationary model, the inflaton must go through the oscillatory 
phase after the end of inflation. As has recently been pointed out, if the gravitational fluctuation 
has inflaton dependent mass term, there will be a resonant amplification of the amplitude of the 
gravitational wave during the oscillatory phase of inflaton. Because of this large enhancement 
of the amplitude of the gravitational wave due to parametric resonance, we show that universe 
can naturally go through the pre-reheated phase with minimally coupled matter field. Therefore, 
during the reheating phase, there is no need to introduce any arbitrary coupling between the matter 
field and the inflaton . 
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I. INTRODUCTION 


fl 


During the last decade, observational cosmology [1[ has reached to a level, through which 
we can go backward in time, and see our universe with a great precession almost near the 
beginning of our universe. By now, we have understood very well the mechanism of how 
the observed very tiny fluctuation {6T/T ~ 10“^) in temperature in the Cosmic Micro-wave 
Background (CMB) sky plays the roll of seed held for the structure of the current universe. 
So far inhation, which is believed to have happened within a very short period of time after 
the big-bang, is believed to be the best theoretical explanation of the origin of those tiny 
huctuation. The inhation was hrst proposed {2 4| to explain some important conceptual 
problems persisted in the standard hot big-bang cosmology. In spite of its simplicity, it was 
soon realised that constructing such a model within the known theory is very difficult. Main 
problems in constructing an inhation model are the mechanism of how it has happened, the 
end of inhation, how the universe gets reheated after such a super cooling phenomena, and 
most importantly the origin of inhaton itself. In addition to all those, there generally exist 
many issues within a particular model under consideration. Living with many fundamental 
questions, in this paper we will ask the question as to how reheating happened after the 
end of inhation. Immediately after the proposal of inhation, this question had naturally 
emerged and studied in detail. There exist lot of studies on this particular mechanism 
{ 5 , 6 ], which we generally call pre-heating. In this paper we will introduce a new mechanism 
in the framework of a certain class of gravity theories where the gravitational huctuation 
acquire inhaton held dependent ehective mass term either from the original theory such as 
massive gravity theory or from some higher derivative gravity theory. In this paper we will 


Q. 


be focussing on a class of massive gravity theory proposed in [7i. In the aforementioned 
class of theories our mechanism will be generic in the sense that for any inhationary model, 
after the end of inhation, the inhaton held generically oscillates, which in turn will provide 
the oscillating mass term for the graviton, and resonant amplihcation in the matter sector 
will naturally happen without invoking any further coupling parameters. In order to realise 
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this, in this paper we will consider a class of ghost free massive gravity model 17|, where 
we only have one extra degree of freedom in addition to the usual two tensor polarization 
modes. And we will identify the additional degree of freedom to be the inhaton. We will see 
that in this class of models the gravitational huctuation gets mass which is dependent on 
that extra scalar degree of freedom. The general physical mechanism to transfer the energy 








during pre-heating period is to couple the matter held as a probe with the background 
oscillating inhaton held. In the usual preheating scenario, one generally couples the matter 
helds with the inhaton held in an ad hoc manner, and try to tune those unknown coupling 
constant to have explosive particle production through parametric resonance. Here we will 
show that because of inhaton dependent mass term for the graviton huctuation how the 
universal minimal coupling of the matter held with the gravity is sufficient to reheat our 
universal after the end of inhation. The idea was hrst put forward in the reference 8| , where 
the inhaton held is assumed to be coupled with the higher derivative gravitational Chern- 
Simons term. However, as we have explained in [^|, it is very difficult to have successful 
preheating with the Chern-Simons coupling. It was shown to be crucially dependent upon the 
large initial amplitude of the gravitational huctuation during reheating. However, this large 
value of initial amplitude is found to be very difficult to achieve. Thanks to the reference 
9|, which pointed out the possibility of generating large amplitude of the gravitational 
wave by parametric resonance considering the inhaton held dependent mass term of the 
graviton. Purpose of the proposal was to create large scalar to tensor ratio without violating 
Lyth bound 10(|. We are interested in preheating. Therefore, to achieve that we will have 
sequence of parametric resonance phenomena. We will see how the parametric resonance 
of the gravitational wave naturally leads to the parametric resonance in the matter sector 
though gravitational coupling. In this paper we will try to consider a minimal model to 
realise this mechanism. 

We structured our paper as follows: in the next section-!! we will introduce a class of 
massive gravity theory, which we will be studying. In the subsequent section-HI, we will do 
the linear perturbation analysis in the gravity and the matter sector. As there in no direct 
coupling of the matter held with the inhaton, we need to take into account the higher order 
term in matter held equation. In section-IV, we numerically solve those equations of motion 
and show that how the non-linear interaction term between the gravitational huctuation and 
the matter held is responsible for the explosive particle production. The main driving force 
behind our mechanism is the successive resonance phenomena. The hrst resonance happens 
in the gravitational sector because of the oscillatory mass term and the resonance in the 
matter held happens in the resonant gravitational background. Because of this resonance in 
a resonant background, the particle production turned out to be much more efficient in our 
model than in the usual re-heating scenario. In the hnal section, we will make some hnal 
remarks on our mechanism and also discuss about some future works. 








II. MASSIVE GRAVITY MODEL 


Massive gravity theory has a long history starting from Fierz-Pauli [m|, who hrst pro¬ 
posed a linearised theory of massive spin two held. Subsequently it has been shown that the 
linear theory has some pathological behaviours which can not be cured in the linear regime. 
Therefore, non-linearity has been introduced at the action level in order to have a consistent 
interacting spin two held theory 12(]. One of the interesting facts of the non-linear modih- 


cation proposed in 


12| is its near unique structure when the Lorentz invariance is imposed. 


But subsequently, the theory has been proved to be difhcult to make phenomenologically vi- 


able because of its low euetgy cut off M. Howevet, utodiffcauttou cau be .uade by suttably 
modifying the Vainshtein scale to extend its regime of validity [1^. Constructing model 
would be easer if one relaxes the Lorenz invariance. Therefore, further modihcation has 
been proposed to violate Lorentz invariance in a controlled manner such that the theory 
remains well dehned and also becomes cosmologically viable j?. 
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16| . As we have already 


mentioned in the introduction, in this paper we will be studying a particular class of models 
proposed in [rj. We want to demonstrate our mechanism considering the minimal model 
where we have gravity with one inflaton held. In order to achieve this one considers the 
following internal symmetries in the goldstone boson sector of gravity 


= = ( 1 ) 

where, Rij is the rotation matrix in the held space, and is arbitrary function of (fP. 

With this assumption we can write most general action for massive gravity as follows 


C = MlR + C{X, 


( 2 ) 


where, £(X, is the most general action for the goldstone boson helds (0°, 0*) in accord 
with the aforementioned internal symmetries. This potential term will contribute to the 
mass term for the graviton in the unitary gauge. Based on the aforementioned residual 
symmetries in the goldstone sector, one can write the following scalar quantities which we 
introduced above in £(X, 
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Since the theory has the symmetry 0'* = 0* -|- F*(0°), by using that one can remove three 
propagating degrees of freedom from the 0® with its hxed vacuum expectation value (0®) = x®. 












By using the Hamiltonian analysis also it has been shown that the above theory has only 
three degrees of freedom. We identify the extra degree of freedom (jp = (p as an inflaton in 
our subsequent analysis. For the sake of our numerical computation we will consider the 
model of inflation from our previous paper 17|]. In compatible with the above symmetries 
we consider the following Lagrangian, 
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(4) 


where, we dehne 

To the leading order, the crucial difference between the above Lagrangian and the usual 
massive gravity theory in the literature is the extra scalar held dependent mass term {q]. 
We are not going to study the theoretical aspects of this kind if mass term in this paper. 
As one can easily identify, this is the term which plays a crucial role for our mechanism 
to work. According to our previous discussion, we will study phenomenological aspects of 
this term in the cosmological context. If one considers the action in unitary gauge, the last 
term in the above action will play like the mass term of the graviton, s is the dimensionful 
parameter. M(0) is the some specihc dimensionless function of the inhaton held. / is the 
decay constant. A can be identihed with the scale of inhation. It is also interesting point to 
note that the mass term will not play the role in the background dynamics. 

We would like to emphasize that in this work we are only interested in the regime where 
the inhaton held coherently oscillates after the end of inhation. To have successful reheating 
this is one of the important requirements for any model of inhation. Therefore, our hnal 
conclusion of this work is valid for any model of inhation. The new ingredient that we have in 
the Lagrangian is the extra term related to the mass of the graviton. The extra mass term is 
constructed in such a way that it only contributes to the mass of the gravitational huctuation. 
Similar situation happens for gravitational Chern-Simons term coupled with the inhaton. 
The crucial diherence is that the Chern-Simon term provides oscillating anti-dissipation 
term for the dynamics of the graviton. On the other hand here we have the oscillating 
mass term. Thanks to the reference , in which the authors have already shown that in the 
oscillating inhaton background the amplitude of the graviton huctuation undergoes resonant 
amplihcation. In this paper will show how that resonant amplihcation helps to have resonant 
amplihcation of the matter held through minimal gravitational coupling. In the subsequent 
analysis we will show how this mechanism works. 






III. TENSOR PERTURBATION AND PRE-HEATING 


In this section we will analyse the tensor and matter perturbation in the oscillating infla- 
ton background. As is well known, in the usual pre-heating scenario, in order to create the 
matter held from the inhaton, one invokes different types of arbitrary coupling parameters 
in the matter Lagrangian, 

Cmat = ]^d^0d^0 - H- (6) 

where, x ^md g are the coupling parameters. In the usual Einstein gravity, since there 
is no parametric resonance in the gravity sector after inhation, one needs to introduce 
above coupling parameters for every matter held to obtain the oscillating mass term for the 
matter held under consideration. Parametric resonance will occur because of those coupling 
parameters. This in turn will provide us the explosive particle production. What is novel 
of our mechanism is that, for the purpose of explosive particle production, we do not need 
to introduce any arbitrary coupling constants for every matter helds. We will show in our 
following analysis, how minimal gravitation coupling is sufficient for this purpose. The 
only parameter that we need in our study is the gravitational mass m(0), which we have 
introduced before in the massive gravity action. Therefore, for our purpose, we set 


1? = 0 ; x = 0 (7) 

As we have already mentioned, at the linear order in perturbation, the extra mass term of 
the gravitational action will only contribute to the mass of the gravitational huctuation. 
The background evolution is governed by the inflaton held 0 with the following metric 

ds^ =—dt^ + a{ff{dx^ + dy^ + dz^). (8) 


After the end of inhation, the inhaton start to oscillate which is known to be responsible for 
the reheating of our universe. For current study we will consider the solution of the oscillating 
inhaton from our reference 17| as shown in the hg{TJ As we have emphasized before, our 


mechanism does not depend on the particular model of inhation under consideration. As 
long as we have oscillating background inhaton, we will have broad parametric resonance 
for the gravitational huctuation for a wide range of momentum modes depending upon the 
oscillation frequency of the inhaton and the coupling paramter A. At the linear order, the 
scalar held huctuation will follow the same dynamics as before. Since our mechanism of 
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FIG. 1: Background evolution of the inflaton field cj) in cosmological time t. We measure the value 
of (p in unit of / and cosmic time t in unit of s. 


pre-heating is mainly dependent on the dynamics of gravitational fluctuation, we consider 
tensor perturbations as follows 

ds^ = —dt^ + a{ty{6ij + hij)dx‘^dx^ (9) 

Apart from inflaton we only have two tensor polarization degrees of freedom. Therefore, 
we can always choose the following transverse and traceless gauge condition for the tensor 
perturbation, 

= 0 ; = 0 . ( 10 ) 


The fourier mode function of hij defined as 


hij{t,x) = 


d^k 




( 11 ) 


Where, efj{k) are the two independent polarization tensor of the gravitational wave. To the 
linear order in cosmological perturbation theory, the Fourier transformed gravitational wave 
equation turns out to be 

k + 3Hh^ + (^ + = 0 - ( 12 ) 

This equation is modified Mathieu equation with the oscillating inflaton background. The 
usual re-heating phenomena is based on this same form of equation of a matter held. Analytic 
study of this equation has been done quite extensively in the context of pre-heating. The 
important difference here is that, resonant production of gravitational wave initiate the 
pre-heating process through the minimal coupling of the matter held. 

By using the properties of the gravitational wave, the form of the scalar held equation 
becomes 


1 

9 + 3He - + —dAe + mle = o. 


(13) 








As one notices that in order to have gravitational effect we restrict np to the second order 
term. Considering the following Fonrier transformation for the scalar field, 


9(t, x) 
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+ h.c.], 


(14) 


the mode eqnation for d(k), take the following from 

i + 3 hI + (F + mjj e = - ^ IY, 


k') d(k') d^k'. (15) 


In the integral part of the above evolution equation for the scalar mode, we consider the 
evolution of a particular scalar mode of k. For a fixed vector k, one can parametrize the 
two tensor polarization modes of the gravitational wave in the above integral as. 


e*/(q) = el(q)4(q)+ 4(q)el(q), 
e 2 ‘'(q) = eUq)ei(q)-4(q)4(q)- 


where, in a particular basis we write 
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(16) 


considering q = (k' — k). As we have mentioned before, for a fixed scalar mode dk with 
momentum k = kz, if we parametrize the q = {k' sin a cos f3, k' sin a sin f3, k — k' cos a}. With 
the above parametrization, the equation of motion for the scalar mode function turned out 
to be 


9 + 3H^ + 
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k'^ sin^ a h(k — k') 
k 2 _|_ /j /2 _ 2kk' cos a 


9{k') dadk'. 


(17) 


Now this is a modified Mathieu equation, where the oscillatory background is set by all 
possible modes of the gravitational wave h(k — k') , which satisfy the following equation , 


+ 3Hh^ + 


k"^ + k'"^ — 2kk' cos a 


+ = 0 . 


(18) 


In order to understand the behaviour of these equations, we will consider two different limit 
on the possible value of the scalar mode momentum. First let us consider large k limit 
compared to k'. The modes in this range will contribute the most. In this limit, the solution 
of h(k — k') becomes approximately independent of k'. Therefore, the above equation for 







the scalar mode function will take the following simplihed form after performing the angular 
integral, 


e + 3H9 + 


+ ml 


fk« 8(1,') ik'. 


(19) 
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To understand the scale around which we will have the resonant amplihcation within the 
perturbative limit, let us introduce the natural time scale coming from the period of 
oscillation of the inflaton held after the end of inhation as mentioned in hg{TJ If we measure 
our time in term of s unit, then the above scalar mode equation turns out to be 
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( 20 ) 
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For the purpose of our analytic discussion, we set the scale factor a(t) = 1. This is a 
reasonable assumption during pre-heating, as the scale factor remains almost constant. We 
also use this fact for our numerical analysis in the later section. At this point we also would 
like to point out that the combination s^h(k) is dimensionless as can bee seen from the 
fourier expansion of the gravitational wave eq. flTT]) . Therefore, our perturbative analysis 
will be under control as long as one satishes s^h(k) < 1. From the usual properties of 
Mathieu equation, the broad parametric resonance will occur if the coefficient of 0(k') in the 
integral part of the above scalar mode equation is larger than unity. Therefore, from the 
scaling properties of momentum k' and h(k), we can set the following qualitative condition 
between the magnitude of the gravitational held and the momentum modes contributing to 
the scalar particle production as follows, 

V\ 5 . „ s 5 
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( 21 ) 


Furthermore, within the perturbative regime, the modes which will contribute to the broad 
parametric resonance for the scalar amplitude should satisfy the following condition. 


k' 


< s^h(k) < 1. 


( 22 ) 


However, if we consider, k' to be larger than k, following the same procedure we have 
discussed so far, the equation of motion for the scalar mode function will turn out to be 
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Following the similar argument as we have discussed so far, the constraint on the momentum 
modes will be as follows 


s3h(k')^(k') 


dk' 


(23) 
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At this point would like to emphasized the fact that the conditions for the momentum modes 
eqs. fl22l) and eq. (l2^ for two different limit of the scalar momentum k » k' and k « k' 
respectively, are within the perturbative limit. The usual parametric resonance will lead the 
gravitational amplitude to be in the non-perturbative regime. In our present analysis we will 
confine ourselves in the perturbative regime by limiting the amplitude of the gravitational 
wave to be less than unity. However, non-linear analysis must be needed to understand 
complete dynamics of the system. Furthermore, it may also be possible that back-reaction 
of the scalar particle production into the gravitational wave dynamics keep it always in the 
perturbative regime. We defer all these issues for our future studies. 

We have not introduced any direct coupling of the matter field with the inflaton. There¬ 
fore, it important to note that in the dynamics of the matter field 9, the gravitational 
fluctuation will play the role at the first non-linear order. However, in order to realise our 
mechanism, not only we need to choose suitable mass term m((/)) for the graviton fluctuation, 
but also we need to be careful about its effect during the inflation. In the usual massive 
gravity theory, modification is tractable in the infra-red limit as the theory remains pertur¬ 
bative in nature. Non-linearity of the theory depends crucially upon the mass term of the 
graviton. In the current construction, we consider gravitational mass is dependent on the 
inflaton field itself. Therefore, one needs to study the strong coupling aspects of the theory 
which we will defer for our future study. Motivated by the reference j9| we assume, in order 
not to have large power suppression of the gravitational wave during inflation, 
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m(0)^ = 


1 + 


0 * 


4 ’ 


(25) 


where 0* is the value of 0 at the end of inflation. Such a functional form will lead to very 
small mass of the gravitational wave during inflation, since inflaton field will have very large 
constant vacuum expectation value during inflation. One can easily compute the tensor 
power spectrum during inflation with almost constant mass limit of the graviton to be 


Ph = 


2H^ 
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(26) 


Therefore, as one can see that due to the specific choice of our mass funtion m(0), modified 
part of the tensor power spectrum due to the mass term is highly suppressed. However, after 
the end of inflation, since the inflaton undergoes damped oscillation around the minimum 
of it’s effective potential, the amplitude of the inflaton will be very small compared to 0*. 
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FIG. 2: Behaviour of gravitational wave Re[hk] in unit of the matter field 0k and particle 
number density re^ in time t which is again measured in unit of s. We have plotted for kx = 
Q,ky = 170,^2 = 130 in unit of s. We chosen the dimensionless graviton mass parameter A = 1. 
We also checked that as we decrease the initial amplitude of the gravitational wave, the resonant 
production of the gravitational wave happens at later time which leads to the particle production 
at late time. 


Therefore, during oscillating phase, we can approximate, 

m{(j)Y = (27) 

where A is dimensionless coupling constant, observing the above expression for the inflaton 
dependent mass term, one would expect to get back the usual massless gravity theory long 
after the end of inflation. For a class of massive gravity theory [funder consideration, it has 
been shown that it reproduces usual gravity theory in the massless graviton limit. In this 
paper we will be restricted to a particular class of minimal massive gravity models where 
the number of degrees of freedom is three. 


Now considering the gravitational wave equation with the above form of the gravitational 
mass term 
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= 0 . 


(28) 


As usual, we consider the time scale in unit of s. As has been extensively discussed in 
and we also have observed from our numerical analysis the successive enhancement 
of the amplitude of the gravitational wave happens at the time of zero crossing of the 
inflaton field, and it oscillates with almost the same amplitude in between the two successive 
crossings. Therefore, analytic treatment of this equation suggested in jo], can be done by 
























the method of successive parabolic scattering. The stochastic resonance certainly will lead 
the perturbative gravitational amplitude into the non-perturbative regime. Therefore, one 
needs to go the non-linear regime. However, in the present paper we will restrict our study 
in the perturbative regime such that the approximate energy density of the gravitational 
wave pI 


ph ^ [ dk (29) 

TT^a^ Jo 

where, (3k is the amplitude of the negative frequency part of the full solution, is smaller 
than the inflaton energy density. Because of resonance, the approximately expression for 
\(3k\‘^ ~ where the exponent pk, dehned for a particular instant of zero crossing of 

the inflaton held is given as 


Pk{t) = —In (l + 2e — 2sin0e 2 ^/l + . (30) 

The expression for and 9 is random phase within (0,27r). Therefore, on can 

see the exponential enhancement of the gravtitional held amplitude. In order to be in the 
perturbative regime, we have to consider ts begin very small. 

One may notice, in our current analysis the eq. fl36p .in the dynamics of the gravity wave, 
we ignored the back-reaction of the matter held. It can be easily understood that at least at 
the initial stage of the pre-heating process, the back-reaction should be small. But it must 
be very important after sufficient number of matter particle produced. The full analysis 
taking into account the back-reaction and rescattering of the particles will be considered in 
our subsequent publication. However, as we pointed out, the matter particle production will 
happen at the non-linear level. Thanks to the parametric resonance that has happened in 
the gravitational sector. Because of this resonance, we will show, how efficient the particle 
production is in onr mechanism compared with the usual reheating mechanism. In order to 
quantify the amount particles that has been produced for 6*twith momentum k we use the 
following well known expression of particle number density |6| 

+ |9(k)N - f (31) 
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Where, ook is the time dependent frequency of a particular mode k. 
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where kphy is the physical momentum. In the above derivation of the number density, we 
have used the following canonical quantization condition for the scalar held mode, 

e*{t)e{t)-e*{t)e{t) = i ( 33 ) 

In the next section will study numerically two equations of motion and try to compute the 
particle number density for a particular mode. For simplicity of our numerical computation, 
we hx the value of ky and do the numerical integration for the mode equation of 6k along 
the kz- 


A. Numerical analysis and particle production 


As we have already discussed before, we are only interested in the oscillatory regime 
of the inhaton, our hnal conclusion does not really depend on a specihc model of in- 
hation. In order to get the background of oscillatory inhaton for our numerical cal¬ 
culation, we consider our specihc model of modihed natural inhation [^. We choose 
M(0) = sin^x [l — cosx sin^ x] where x = (p/f. For this specihc form of the function, we 
had the number of e-folding J\f = 50, and the value of the spectral index is Ug = 0.96. In order 
to achieve the aforementioned cosmological quantities compatible with the observations, we 
found the values of the other parameters to be / = 0.84Mp, s = 2.25 x 10“®Mp, A = O.OllMp. 
For our numerical purpose we set f = 0 as the beginning of the inhation. After the end of 
inhation at around f ~ 4 cosmic time in unit of s, the inhaton held will start to have co¬ 
herent damped oscillation as shown in the hg-[Tl At this point, we are again emphasizing 
the fact that our analysis and the hnal result is conhned within the oscillatory regime of the 
inhaton held. The oscillation after inhation is generic to any inhationary model. Therefore, 
our qualitative results are insensitive to the specihc model of inhation under consideration. 
For the sake of our numerical analysis, at time f = 6 we normalize the cosmological scale 
factor a(6) = 1. For a wide range of initial amplitude of the gravitational wave, we found 
the resonant production of the gravitational wave. More over this resonant production of 
gravitational wave will trigger the pre-heating phase of the universe through the matter 
helds with the universal gravitational coupling. In the inhationary background, to quantize 
a held we choose the usual Bunch-Davis vacuum. 


d(k) = lim 
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FIG. 3: Comparing the particle production for matter field 6{k) in massive gravity scenario (left 
panel) with the usual pre-heating scenario (right panel) considering Cint ~ 9 coupling term. 

Plot for the particle density in time t which is measured in unit of s. We take kx = 0,ky = 
170, kz = 130 in unit of s. Choose X = g = l,mg = 0.1s. 

In the above expression, we also assume approximately the de-Sitter vacuum. As we see, 
lower the value of a particular mode k, larger the initial amplitude in the quantum vacuum, 
leads to the earlier production of matter particle of that particular modes after the end of 
inflation. As a particular mode starts evolving from its quantum vacuum dehned in the far 
past and deep inside the Hubble volume during inflation, its amplitude will start to decrease 
due to Hubble friction and hnally gets frozen after the horizon exit. Therefore, during 
reheating the amplitude at the horizon exit will set the initial condition for a particular 
mode of the gravitational wave to evolve. From the aforementioned arguments, we see that 
the parametric resonance for the matter held will happen mostly for the long wavelength 
modes compared to the Hubble length. We can clearly see from the hg-12] that because of 
the resonant gravitational wave production triggered by its oscillatory mass term, minimal 
coupling of the matter held with the gravity is sufficient to reheat the universe. 


we also want to compare our result with respect to the usual pre-heating scenario jb]. At the 
linear order in perturbation theory, in our current scenario we want to see how efficient the 
particle production is in comparison with the usual pre-heating scenario jb]. One can clearly 
see the diherence between the distinct ehect of parametric resonance on the gravitational 
huctuation h{k) and the matter held huctuation 9{k) from the hgj2j Important point to 
mention that in the scalar held sector, the parametric resonance happens in the resonant 

























gravitational field background. Therefore, one can clearly see from the left panel of Figj3] 
that the produced density of particle in the massive gravity framework is behaving like 
double exponential function in time. The initial condition, we have used in our analysis is 
that there in no particle at f = 6 cosmic time in unit of s after the inflation. After the end of 
cosmic time t = 7.64, so that the amplitude of the gravitational wave is in the perturbative 
limit, one can see from the £gj3]that the number of particle produced per unit volume in 
the massive gravity scenario is ~ times larger than that of the usual pre-heating scenario 
for a particular mode. Therefore, in the massive gravity scenario, the particle production 
is much more efficient compared to the usual pre-heating scenario. It is possible that most 
of the energy of the inflaton field may easily be transferred through successive parametric 
resonance hrst to the gravitational sector and then to the matter sector within a very short 
period of time, and the usual reheating follows. The whole mechanism of particle production 
that we have discussed so far is highly non-equilibrium in nature. Therefore, the elementary 
theory of reheating is still necessary to thermalization all the produced particle, and the 
universe to enter into the usual radiation phase. Now, because of highly efficient particle 
production compared to the usual pre-heating scenario, if most of the inflaton energy is 
transferred to the relativistic matter particle, the usual standard model couplings among 
the produced particles may be sufficient to thermilize them. Therefore, we may not need to 
introduce any arbitrary coupling at all to have the successful reheating of our universe after 
the inflation. Of course our current analysis is in the linear regime, and we have not taken 
into the back-reaction and the re-scattering of the fields. We differ the detail studies of all 
the aforementioned points for our future work. 


IV. CONCLUSION 

In this paper, we have tried to constructed a phenomenological model of pre-heating by 
introducing inflaton held dependent mass term for a certain class of massive gravity theory. 
In that particular class of theories, we have three degrees of freedom because of some special 
internal symmetries. For the purpose of our study, we identihed the extra degree of freedom 
to be the inhaton held. During inhation, gravitational mass term will suppress the power 
in the gravitational wave spectrum. In order to avoid large suppression, we parametrize the 
mass function in such a way that during inhation, its value is very small. On the other 
hand, after the end of inhation, the amplitude of the gravitational wave will be amplihed 


due to parametric resonance because of the mass term[9j. Now the question naturally arises 
is why have not we observed this large amplitude gravitational wave yet? Natural answer 
would be that this large amplitude gravitational wave will contribute largely to reheat our 
universe after the inflation. In this paper we have shown that how this resonant gravitational 
wave background can naturally lead to the pre-heating phase of our universe. Interestingly, 
to pre-heat the universe we do not need to introduce any arbitrary couplings between the 
visible sector matter held and the inhaton held in an ad hoc manner. The natural universal 
gravitation coupling is sufficient to reheat the universe. At the linear order in perturbation 
theory, we also compared the efficiency of the particle production between our mechanism 
and the usual reheating mechanism j^. We found that because of the resonant gravitational 
background, in our mechanism the particle production is highly efficient compared to the 
usual scenario. This phenomena of resonance in a resonant background has not been studied 
before at least in the cosmological context. It would very much interesting to study this 
mechanism in detail in various physical situations. 

Using the same mechanism that we discussed so far, we can study the production of chiral 


gravitational wave c 


derivative term 


18|, 


uring preheating if we further introduce the Chern-Simons (CS) higher 
19| coupled with the inflaton held in the action, i.e. 


Ccs — j4>RR‘ 


(35) 


The inhaton held is now like a pseudo scalar axion, where 


RR = R\^f^ R^ 


flap 




i/7<5- 


is the usual 4-dimensional Levi-Civita tensor, g is the free coupling parameter. Since in 
the background of inhaton this term violates parity, the left and right handed gravitational 
wave propagate diherently according to the following equation 


Ml + 

^ fa 


+ mh\ + -h\] +m{ct>yh\ 


k g va 
f a 


(^j, - H4,) , 


(36) 


where A = L,R are left and right circular polarization mode of the gravitational wave. 
Therefore, we will have resonant chiral gravity wave production. This ehect may be very 
important to produce the leptogenesis in the early universe. As is well know 2^, if we 


couple gravity with the leptons, the total lepton number current will no longer be conserved 











because of gravitational anomaly, 


d,r = 


IbTT^ 


< RR >, 


(37) 


where, < RR > is the quantum expectation value of the Chern-Simons term. Therefore, 
non-zero expectation value of this Chern-Simons operator will produce the leptogenesis. 


This mechanism has already been discussed in 


18| . However, in our model the resonant 


production of chiral gravitational wave will happen during oscillation period of the inflaton. 
Therefore, the maximum lepton number violating contribution will be coming during the 
preheating period. Currently we are looking into this issue in greater detail. We are also 
looking into the possibility this new reheating mechanism in a Lorentz invariant higher 
derivative gravity theory coupled with inflaton held. We know that the higher derivative 
term such as Gauss-Bonnet term coupled with the inhaton held contributes an ehective 
inhaton held dependent mass term in the gravitational huctuation. Currently we are also 
looking into that possibility. 

We have also pointed out that the specihc class of massive gravity theories, we have 
considered for our study, need to be studied in detail. In massive gravity theory, one generally 
introduces the constant mass term for the graviton. But here we have introduced extra held 
dependent mass term 9j. This particular held dependent mass term should be studied with 
a special care. As one knows interacting massive spin two held coupled with other held such 
as our scalar held has several constraints to be satished. We keep all these questions for our 
future studies. 
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